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Abstract 

> . 

\ We present a stochastic theory of charges moving in an electromagnetic field us- 

ing nonequilibrium quantum field theory. We give a first principles' derivation of the 
Abraham-Lorentz-Dirac-Langevin equation which depicts the quantum expectation value 
for a particle's trajectory and its stochastic fluctuations by combining the worldline path 
integral quantization with the Feynman- Vernon influence functional or closed-time-path 
effective action methods 0, §]. At lowest order, the equations of motion are approximated 
by a stochastic Lorentz-Dirac equation. 
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1 Introduction 

1.1 Particle versus field formulation, fixed versus dynamic back- 
ground 

The contrasting paradigms of particles versus fields give very different representations for phys- 
ical processes. The success of quantum field theory for the past century has propelled the field 
concept to the forefront, with the particle interpreted as excitations of field degrees of freedom. 
In recent years, the use of the quantum-mechanical path integral in string theory has inspired 
a renewed effort towards particle-centric quantum formalisms. The so-called worldline quanti- 
zation method - where the particle's spacetime coordinate x M (r) is effectively quantized - has 
been used for calculating high-loop processes in non-dynamical classical background fields ||. 
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The background field approximation is useful and well-defined in conditions where quantum 
fluctuations of the fields are small. However, many problems require the inclusion of field dy- 
namics beyond the fixed-background approximation. A simple example is radiation reaction. 
More examples are found in semiclassical gravity, where spacetime playing the role of 'particle' 
is treated classically in its interaction with quantum fields |4[]. 

The proper treatment of radiation- reaction requires, at the least, consistency between the 
particle's (quantum) average trajectory and the average radiation. When one treats the particle 
as quantum mechanical, but coupled only to the mean-field background (and not the quantum 
field fluctuations) one is working within the semiclassical regime ||. (When quantum fluctua- 
tions of the particle motion are sufficiently small, one recovers the classical regime 0). In our 
program, we go beyond the semiclassical approximation to include the influence of quantum field 
fluctuations on the moving particle J7|, thus entering the stochastic regime. To do this properly, 
we must begin with the full quantum theory for both particle and field, since a particle can 
never fully decouple from the field. Also since most field degrees of freedom remain unobserved, 
we treat the particle as a quantum open system with the field acting as an environment J8|. 
This is how nonequilibrium and stochastic mechanics ideas enter, and why the self-consistent 
quantum particle evolution is generally described by stochastic quantum dynamics 0. 

The open system concept and the accompanying method of coarse-grained effective action 
and influence functional have been applied to interacting quantum fields and semiclassical 
gravity. The new task we undertake is the adoption of the worldline quantization method for 
relativistic charged particles treated as an open system [0, [K| . This is an alternative formulation 



of QED in terms of quantized worldlines and quantum fields. The worldline structure is highly 
efficient for describing the particle-like degrees of freedom (such as spacetime position), while 
the field structure most easily describes processes like radiation. 

1.2 Radiation reaction, causality, and the Abraham-Lorentz-Dirac 
equation 

The conceptual and technical scheme we construct is advantageous for carrying out calculations 
which conventional methods find difficult, and for clarifying confusions in existing theories. 
Consider the classical problem of radiation reaction. A relativistic solution to classical radiation 
reaction is known as the Abraham-Lorentz-Dirac (ALD) equation, whose solutions give rise to 
famous paradoxes including: runaway solutions, pre-acceleration, and the need for higher- 



derivative initial data [11]. These difficulties with the classical theory may be avoided by 



replacing point particles with extended objects p!2 |; but a relativistic treatment is cumbersome 



Moniz and Sharp |L2] have shown that non-relativistic QED for extended objects can be causal 
and runaway free. Low has shown that runaways apparently don't occur in spin 1/2 QED 



[ T3| . Others |14j] have shown how the lowest order-perturbation expansion of the ALD equation 
can be inferred within conventional relativistic-QED field theory. But most of these earlier 
works lack a comprehensive theoretical foundation for addressing all these issues of quantum 
field theoretic and nonequilibrium (in fact non-Markovian) nature, and for dealing with their 
interrelations consistently and completely. We show, for instance, that the equations of motion 
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for the particle's mean trajectory are causal and runaway-free regardless of fine structures like 
spin. We highlight the role of the field-environment in decohering the particle trajectory. Going 
beyond the result in [14| which is to order e 4 in the particle charge, we find the ALD equation to 
all orders in e for the semiclassical limit. This shows the significant advantage of the worldline 
formulation together with the loop expansion in the influence action. The inclusion of spin and 
color are important to make full use of this approach in QED and QCD. 

Beyond demonstrating the conceptual and technical consistency of the particle equations 
of motion, we describe the influence of quantum field fluctuations on the particle trajectory. 
Quantum field fluctuations play dual roles: they are responsible for the decoherence in the 
quantum particle (system), leading to the emergent semiclassical and classical behaviors, and 
they provide an effective classical stochastic noise in the particle motion |J . We show that the 
one loop equations of motion for the particle trajectory are approximated by stochastic differen- 
tial equations featuring colored noise that encodes the influence of the quantum statistics of the 
field . When there is sufficient decoherence, these equations give excellent approximations to 
the particle motion. Even in the case of weak decoherence they provide a one-loop approxima- 
tion to the particle correlation functions, and may be extended self-consistently to higher-loop 
approximations via the N-particle-irreducible (NPI) effective action [HJ. In the second part of 
this paper, we show the main steps in our derivation of the ALD-Langevin equations describing 
the stochastic dynamics of particle motion. 



2 Coarse-grained effective action and equations of mo- 
tion 

2.1 The in- in generating functional 

We begin with the construction of a generating functional for the particle worldline correlation 
functions. We assume an initially factorized density matrix for the particles plus field at initial 
time ti with the form 

P (ti) = Pparticles <g> p fieM (ti) = f) z (U) <g> f>^ (t { ) . (1) 

For simplicity, consider only one particle worldline and assume that the initial particle state is 
a positive frequency configuration space state p z (U) = \z\ )(z\ |, where \z\ ) = +) = 

(zj) \vac) particle- These are quasi-localized, physical, relativistic one-particle states; we treat 
more general initial states in |7|, [Hj. A basis of field states is \A (x)); we define a direct product 
basis for the particles plus field by \zA) = \z) <g> \A). 

To find the generating functional for the correlation functions defined as 

(n) (r„)) = Tr Az (n) ...z v (r n ) p in ) , (2) 

where z^ (r) are Heisenberg operators and the trace is over all final particle and field states, we 
need the time evolution of the density matrix: 

p(t f ) = U(t f ,t t )p(t t )U(t f ,uy . (3) 
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In field theory, one could introduce a functional representation for the unitary evolution oper- 
ators U in terms of the sum over all field-histories between tj and tf, consistent with the initial 
and final state boundary conditions. 

Instead, in a hybrid particle-field theory we start with the action 

S[z,A] = S g [z]+S A [A]+S int [z,A] 

= J dr (z^z^/N - Nm 2 + (r) z" (r)) - J d 4 x F 2 /A 

d A xf(x;z){A ll (x)+A lx (x)). (4) 



The quantum degrees of freedom are the field A^ (x) , the particle coordinates (r) , and the 
'lapse function' N (r) . The action S [z, A] is invariant under reparametrizations r — ► r' = 
t + e (r) ; this is a crucial symmetry of the relativistic worldline formulation. These issue are 
discussed in detail in J7|. The F 2 /4 term is the usual free electromagnetic field action. We allow 
for the possibility of an additional, non-dynamical background field A. The particle couples to 
the field through the current (x; z) , which is given by 

f (x; z) = e J dTz» (r) 8{x-z{r)). (5) 

In ([D we add a source term (r) whose role will be made clear below. 

The in-in or closed-time-path (CTP) generating functional is given by flC 

Z in _ in [h, ti] = Ti Az (lJ h (t f , ti) piU h < (t f , ti)A 

= J dzfdz'fdAfdAidA'J (z f - z'f) (d' tf - d tf 

x(Af Z f\A4 +) ) h (AfZ , f\A l z < i +) r h/ p(A l ,A l ;t l ) . (6) 

tf=t' f 

We have implicitly assumed in writing dfj) that the sole particle worldline begins at Z{ on the 
initial time hypersurface £ (ti) and ends somewhere on the final time hypersurface S (tf) . We 
then integrate over all possible Zf on E (tf) . A more general form, with more general boundary 
conditions, is given in 0, [TIJ. For example, a worldline that begins and ends on S (ti) (E (tf)) 
describes pair-annihilation (creation). 

We may now use the worldline formalism to give a path integral representation for the 
propagators (AfZf\Aiz\ )h an d (A'^i \Afz'f)h' = (Afz'^A^z^)^,. They are given by 

(AfZf\A l z i ^ ) ) h = [ Z,Ai DzDAex V Us[z,A,h]\, (7) 



n 



and similarly for (A^z- \AfZf)h>- The measure is defined by 

rZjA$ ntj poo p 

i ' DzDA = YiX\. \ dz i / dZi / dA »( x j) 

J ZiA; ■ ,, • Jti J — oo J 
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The generating functional then takes the form 



^-in [h, ti) = J dz f f£f-t DzDz' (d' tf -d tf )e 



,i(S z [z,h]-S z [z',h']) 



"■f= v f 



,AtA' t 



X 



x / ' / dAfdAidA'iDA DA'p (A h A'-, U) 

% i 

exp {i {3 A [A] + S mt [z, A] - S A [A'} + S int [z> , A'})}. 



(9) 



The sum over histories in (0) include all worldlines that begin at Zi, end at Zf, and are bounded 
between the initial and final times, ti and tf, respectively, including those that are spacelike 
and that change direction in time @. 



2.2 Fey nman- Vernon influence functional 

Evaluating the generating functional for our action (f|) is virtually impossible. The interaction 
term Si nt , while linear in A, is highly non-linear in the worldline coordinate z. However, if the 
initial state of the field p v (A iy A'^ ti) is Gaussian, the entire integrand in the second line of (0) 
is Gaussian in the variables A. Thus, we can then do the A, A' path integrals exactly. Initial 
Gaussian field states include thermal, squeezed, and coherent states, and therefore provides a 
fairly rich set of interesting and physical examples. The result of the A, A' functional integrals 
is the Feynman- Vernon influence functional M: 

r A f A 'f 

F[z,z'] = / dA f dA i dA , i DADA'p{A h A' i ;t i ) (10) 

J A f A' 



x exp {S A [A] + S mt [z,A] - S A [A'} + S mt [z', . 



Notice that initial factorizability (see (|l|)) makes this result independent of the initial particle 
state p z (ti) . The generating functional is now 



z in . hl [h, h'\ = fdz f r r zf DzDz' (d' tf - d t ^ 

J J Zj J Z- 



x exp \ % - (S z [z, h] - S z [z', h'} + S IF [z, z']) 



t f =f f 



where the influence action Sif is defined via F [z, z'\ = exp (ASjp [z, z'\) . For the Gaussian 
case, the influence action is known; it is 

Sif = j^j J dxdx'g^j~ (x)G R (x,x')j+ (x 1 ) (12) 
+ l - f [ dxdx'g^j~(x)G H (x,x')j-(x'). 



h 
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The are defined by 

3$ ( x 'i z i z> ) = J dr (V (x-z)± z'^5 (x - z')) , (13) 
Gr (x, x') is the retarded Green's function 

G R (x, x') = Tr A ( [i M (x) , A v (x')j Pa ft)) ft f ) , (14) 
and Gh (x,x') is the Hadamard function 

G H (x, x') = Tr A (x) , A, {x') } PA ft)) • (15) 

The coarse-grained effective action is defined by [0] 

Scgea [z ± , h±} = S z [z, h] - S z [z', ti) + S IF [z*] , (16) 
where h$ (r) = (h, (r) ± fcj, (r)) . 

2.3 Decoherence and the (modified) loop expansion 

To evaluate the non-linear z, z' path integrals, we need to use an approximation method. In 
perturbative scattering theory, one expands e lSlF ' h = 1 + iSiF/h + O (e 4 ) . (Lowest order in- 
teractions are O (e 2 ) because we are working in a in-in formalism.) But, to evaluate quantum 
corrections to the semiclassical result, the loop expansion (or even better, the NPI effective 
action jlSl ) is the most suitable. Observe that 

It is not difficult to show that Ihi(SV,f [z^]) > for all z ± . For large z~ = z — z', it is often 
the case that \F\ <C 1. This suppression of the generating functional for large z~ is a manifes- 
tation of decoherence; it occurs due to the large number of degrees of freedom A (x) that are 
integrated out (traced over) in finding the influence action. We now assume that decoherence 
effectively suppresses large z~ (to what extent this is true depends on the details of the sys- 
tem, environment, and their mutual coupling) and expand Scgea in powers of z~ . We define 
Scgea = ^(Scgea) and S CGEA = Iva(ScGEA)- Dropping the O (z~) 3 terms gives a Gaussian 
approximation to the z~ path integrals. Then the Dz~ integration may be performed giving 

Z in , m [h,h'} = j 'dz f *' Dz + J*'^ Dz-fy-dt^e'S^^W (18) 

x ei 1 ( ££ Sf A ) e ~ I drdT' g ^z-(T)G H (z + (r),z + (r'))z-(r') 

x exp |l J drh-z + » (r) - J drdr\ (r) G H (z+, z +! ) rf (r')| , 
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where we have defined the variable tj^ (t) by 



i^(t) 



SScGEA [ z+ 1 z+ ] 



SZ: 



(19) 



z~=0 



We now have only the z + integration left. Note that ([T^) is an equation for z + as a func- 
tional of 77 (r) . We use this fact to change integration variables from z + — > rj, yielding Dz + = 
(Jacobian) xDr] = Drj, since the Jacobian equals 1. This gives us our final result, the generating 
function 



with 



Z in . in [h~] = J D V P [rj\ exp | % - J dr h~ (r) (r) 



P [V] = exp |- j drdr'r], (r) G H (z + , z + ') ^ (r')\ . 



(20) 



(21) 



P [rj] may now be treated as a probability distribution for rj (r) , thus rendering rj' 1 (r) a stochastic 
noise variable. ( |T9"D is therefore a Langevin equation for z+ (r) with noise 77^ (r) . Derivatives 
5 n -^in-m I ' {t\) ...8h~ {j n ) give the correlation functions @j. Because z + is a nonlinear 
function of 77, as determined by (|T9|), and G// (^ + , ^ +/ ) is a function of z + , the distribution 
P [77] is not Gaussian. In [Kj, we show how to perform a cumulant expansion on the Langevin 
equation. The order e 2 term then generates Gaussian noise, but more generally, one recognizes 
that even after decoherence has suppressed large quantum fluctuations in z~, the effective noise 
will be non-Gaussian for nonlinear Langevin equations. Only when the Langevin equation is 
linear is the noise strictly Gaussian; at high temperatures it is approximately white. 



3 The ALD-Langevin equation 

To apply this general result to charged particles in the quantum electromagnetic field requires 
regularizing Scgea by introducing a high-frequency cutoff A in the Green's functions G\ H . We 
perform this analysis in [[K]]. The resulting Langevin equations of motion found from (|19|) are 



m (r) z^ (r) = ez^ 



+ e 2 X> (n) (r)< ) +?/V('r)- (22) 



F^y is included as a possible classical, external background field. The variable r = r — Tj is the 



elapsed proper-time from the initial time Tj. 

Quantum field fluctuations dress the particle state giving it an effective mass m (r) . At r — 
(corresponding to the initial time when the field and particle are uncorrelated) m (0) = mo, the 
bare mass of the particle. The mass is then renormalized on the cutoff timescale 1/A, quickly 
reaching its late-time value m (00) = m + a A. The constant a depends on the details of the 
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cutoff since the high-energy theory ultimately determines the precise relationship between bare 
and effective mass. 

Radiation reaction is given by the second term on the RHS of (p2]). The 



u f) = U W (i, 2,..., d n+1 z/dT n ) (23) 

are combinations of higher time-derivatives of z(r). The time dependent coefficients (r) 
determined by the regulated retarded Greens function scale as A 2 ~ n with the cutoff. There- 
fore, the n > 2 terms are high-energy corrections to the radiation reaction force and as such 
they are strongly suppressed at low-energies. The late-time behavior of the n = 2 term gives 

/r- (2) = £ (^ 2 + *j > (24) 



which is just the usual ALD results. In the expansion (22) is shown to be convergent 
for all r, including r — ► oo and r — > 0. The short-time behavior is particularly significant: 
lim r ^ 9^ ( r ) = for all n; therefore, the radiation reaction force identically vanishes at £j, and 
the usual Newtonian initial data (i.e. position and velocity) uniquely determine the equations 



of motion. In flOU , it is further shown that the equations of motion are strictly causal and 
without runaways. The non-Markovian time- dependence of the radiation reaction force reflects 
the fact that in any nonequilibrium quantum setting it takes time for the particle's self-field to 
adjust to changes in its motion; it can not respond faster than the effective high-energy cutoff 
A. This lag is always enough to preserve causality. To the extent that QED is a good effective 
field theory, these results demonstrate the consistency of the low-energy physics, including the 
ALD equation as the proper low-energy limit for radiation reaction. 
Finally, the noise is given by 

{{V»{r) ,vAr)})s = 9»vGi{T,T') , (25) 

where () s is the noise average with respect to the distribution P [rj] . Consistency with the 
fluctuation-dissipation relation derived in |L(J requires mutually consistent regularization of 



both Gr and Gh- Therefore, the noise will also depend on the cutoff A. In fact, the noise 
kernel is modified at early times by the restriction that the path integrals over z ± are bounded 
between the surfaces £(£«) and £ (t/) . This modification of the noise is balanced by the early 
time radiation reaction behavior. At late times, the restriction on the paths has little effect on 
the integrals and the noise quickly approaches the form ( pop just as the radiation reaction force 
approaches its late-time limit. 
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